HIGH -ALTITUDE ROBIN DATA-REDUCTION PROGRAM 

By James Luers and Nicholas A. Engler 

University of Dayton Research Institute 
Dayton, Ohio 

SUMMARY 

The problem of computing winds and thermodynamic data utilizing the space -time 
coordinates of a falling sphere becomes complex when the apogeee of the sphere is over 
100 km. This paper describes the methodology used in constructing the computer 

program. 


ROBIN/ ARCAS SYSTEM 


The ROBIN /Areas system consists of a ROBIN balloon, an Areas rocket motor, and 
an AN/FPS-16 tracking radar. The ROBIN sphere is made of J - mil Mylar inflatable to 
a diameter of 1 meter containing an internally supported corner reflector. Packaged in 
a collapsed condition within the nose-cone of a meteorological rocket, it is ejected at the 
apogee of the rocket and inflated to a super pressure of approximately 10 millibars by 
vaporization of a liquid such as isopentane. Thus inflated, the ROBIN sphere is tracked 
from apogee to approximately 30 km altitude by an AN/FPS-16 high-precision tracking 
radar. The Areas rocket motor is a 4.5-inch-diameter solid-propellant end-burning 
rocket capable of carrying the sphere payload to an altitude of 75 km. The FPS-16 
tracking radar generates spherical space -time coordinates at digitized increments of 
1/10 second. From the space-time coordinates, the meteorological parameters of den- 
sity, wind, temperature, and pressure are deduced. A discussion of the ROBIN/Arcas 
system, with some of its advantages and shortcomings, is contained in reference 1. 

Early results from the ROBIN/Arcas system whetted man’s appetite to extend 
the passive sphere experiment to altitudes beyond the reach of the Areas motor. To 
achieve this dream the Air Force Cambridge Research Laboratories (AFCRL) has 
experimented with a variety of boosted rocket motors. The most successful of these is 
the Viper-Dart rocket motor. The Viper-Dart rocket is capable of carrying the ROBIN 
payload to an apogee of 125 to 140 km. It was anticipated that the extended balloon apogee 
of 125 km would enable density and perhaps wind measurements to be extended to 100 km. 

The data-reduction program designed by Engler (ref. 1) to reduce the data from 
ROBIN/Arcas flights produced accurate density and wind measurements below an altitude 
of 70 km. The high-altitude ROBIN/Viper-Dart system, however, produces balloon 
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velocities and accelerations much larger than the ROBIN/Arcas system. For this reason 
the smoothing techniques used in the 1965 ROBIN/Arcas program are not optimum for use 
with the ROBIN/Viper-Dart system. It has been shown by Engler (ref. 2) that the stand- 
ard ROBIN/Arcas data-reduction program is not satisfactory for use with high-altitude 
flights. Accepting the recommendations of Mr. Engler, AFCRL has requested the 
University of Dayton Research Institute (UDRI) to develop a new ROBIN data-reduction 
program which would result in optimum density and wind measurement for high-altitude 
rocket launches as well as the standard Areas rocket. It is the purpose of this paper to 
discuss the new data-reduction program, to explain the rationale and methodology used to 
design the program, and to discuss the errors in the winds and thermodynamic data that 
result from the use of this program. 

PROGRAM SPECIFICATIONS 

The preliminary specifications for the program consisted of the following items: 

(a) the program should be optimum for measuring density and wind in the 70 to 100 km 
region of the atmosphere assuming a balloon apogee of 125 km, (b) the program should 
also give accurate and reliable density measurements from 30 to 70 km, (c) even though 
the data-reduction technique need not be optimum for balloon apogees other than 125 km, 
other balloon apogees between 75 and 140 km should not result in a serious degradation 
of the meteorological parameters, (d) temperature and pressure accuracies should be 
commensurate with density accuracy, and (e) the program should accurately determine 
the altitude of balloon collapse so that density calculations can be terminated. 

DENSITY AND WIND MEASUREMENTS 

To obtain density, the drag force that the atmosphere exerts upon the sphere must be 
measured. In the altitude region from 70 to 100 km, the vertical velocities and acceler- 
ations are much larger than the horizontal velocities and accelerations. For this reason 
the drag acceleration is primarily in the vertical direction. Accurate density calculations 
are thus largely a result of the accuracy to which vertical velocities and accelerations can 
be measured. 

Horizontal winds influence the sphere's trajectory by inducing horizontal excursions 
in its path in three-dimensional space. These horizontal excursions are used to recon- 
struct the wind profile. Thus, for measuring wind, the horizontal velocity and accelera- 
tion components must be determined accurately. 
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Since densities depend primarily on vertical measurements and winds depend pri- 
marily on horizontal measurements, it is possible to construct a data-reduction scheme 
which will optimize both wind and density measurements. The data-reduction scheme 
used to reduce Viper -Dart flights is so designed. 

DENSITY 

Density is computed by the following equation: 

m(g - z - C z ) 

P ±C D Av(z - W z ) + V B g 

The symbols are defined in the appendix. The computed density error is a result of the 
errors present in the parameters on the right side of equation (1). A negligible contribu- 
tion to the error in density is made by Vb, C z , g, m, and A. The remaining vari- 
ables which make a significant contribution to density error are Cd, z, z, W z , and v. 

DENSITY ERROR EQUATIONS 

For the purpose of deriving an error equation for density, the density equation 
(eq. (1)) can be simplified to 

m(g - z) ( 2 ) 

P_ I Cd a(z-W z )2 

where v has been set equal to (z - W z ) and both buoyancy and Coriolis force have been 
neglected. 

Considering the error in density to be a function of the errors in Cd, W z , z, and 
z only and assuming further that the errors in the radar coordinates are independent and 
normally distributed with mean zero and variance a, the error equation for the percent 


error in density is given by 



The object of the computer program as mentioned earlier is to minimize equation (3) in 
the altitude region from 70 to 100 km. Equation (3) cannot be minimized by minimizing 
each of the terms on the right-hand side of the equation because the last three terms are 
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interrelated. The first two terms, however, are independent and thus allow for individual 
minimization. 


ERROR IN DRAG COEFFICIENT 

There is no commonly accepted drag table in use today with known accuracy. Of the 
drag tables being used, disagreements as high as 15% to 20% exist in certain areas (ref. 3). 
As seen in figure 1, disagreement exists not only between the values of Cp for various 
drag tables (refs. 4 to 6) but also in the slope of the curves. The drag tables used in pre- 
vious ROBIN programs designed by UDRI rely primarily on the work of Dr. Helmut 
Heinrich and others (ref. 5). The accuracy of this table is uncertain especially in areas 
where interpolation of the drag coefficient is necessary such as from Mach 0.7 to Mach 1.0. 
In an effort to evaluate the most recent drag table to appear in the literature (ref. 4), UDRI 
reinvestigated the drag-coefficient values of references 5 and 4. The table which appears 
as figure 2 is basically the work of reference 5 in the supersonic region and of reference 4 
in the subsonic region. The impressive aspect of this drag table is the similar shapes of 
the Cj) curves given as a function of Mach and Reynolds number even though the drag table 
was the result of two independent researchers using two different techniques for calcu- 
lating drag. However, even though this drag table shows smooth consistent drag curves, 
it is impossible to quote specific accuracies of the drag table because of the conflicting 
results obtained by the other experimenters and because of the interpolated section of the 
table. The stated accuracies by the experimenters are as follows: 

Goin and Lawrence: approximately 2% 

Heinrich et al. (supersonic): maximum possible error ranges from ±2.3% to 
±27.9%; however, actual errors are usually not the possible maximum 

Since the accuracy of the drag table cannot be determined precisely, it is impossible 
to give an exact RMS error value for the percent error in density when reducing a high- 
altitude balloon flight. It is, however, possible to determine the error in density that 
results from the other terms of equation (3). Improvements and verifications of drag 
results will enable one at some future time to accurately state the true percent error in 
density. 


VERTICAL WIND ERROR 


The density error variance resulting from vertical winds is given by the expres- 
2 


sion 


2oi 


W„ 


z - W„ 


To a falling sphere, a vertical wind looks identical to a change in density. 
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As a result, a data-reduction program cannot distinguish density perturbations from 
vertical -wind oscillations. In order to compute densities, an assumption must be made 
either concerning vertical winds or concerning density perturbations in the atmosphere. 

Assumption A: W z = 0 

Assuming no vertical motions in the atmosphere, equation (1) can be solved by 
substituting W z = 0 on the right-hand side of equation (1) and evaluating all other terms 
by conventional means. Under this assumption, any vertical winds present in the atmo- 
sphere will appear as density oscillations. The relationship between vertical winds and 
density perturbations is exhibited in figure 3 for an escape altitude of 125 km. Care must 
be used in interpreting figure 3, however, because of the smoothing effect produced in the 
program. An example will clarify this point. If a sinusoidal vertical wind varying with 
altitude with amplitude of 5 m/sec is present at 60 km, then this vertical wind would be 
damped by the smoothing and appear in the printout as something smaller, approximately 
2 m/sec amplitude. Since the program attributed the vertical motion to density pertur- 
bations, the result of an actual 5 m/sec vertical wind would, using figure 3, appear as a 
2.4% density perturbation. To effectively determine what vertical wind could have caused 
a density perturbation in reduced data in addition to figure 3 one must know the reduction 
in magnitude of the vertical wind resulting from the smoothing technique applied, i.e., the 
frequency response of the program's smoothing filter to a sinusoidal vertical-wind 
oscillation. 


Assumption B: p = PQe az 

If density is assumed to follow some mean path then perturbations from this path 
can be attributed to vertical winds. Since density varies exponentially with altitude, a 
mean exponential path is appropriate. Using this assumption, vertical winds can be com- 
puted by the equation 


■,, 7 ^th ^emp 
Wz = 2Kz Po e- 


( 4 ) 


A description of the variables in this equation and its application is given in reference 1. 


Since, to the best of our knowledge, meteorologists accept density perturbations at 
least as much as they accept vertical winds, assumption A has been and will be used in 
this program. 
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NOISE AND BIAS ERROR TERMS 


The remaining terms in the density error equation are given by 


2(Xg 

2 

r % i 

2 

z2(z + Az - g) - (z - Az) 2 (z - g) 

Z - w z 

+ 

L z - 6. 

+ 

(z - g)(z + Az) 2 


The first two terms of equation (5) are the result of the noise present in the radar coordi- 
nates (noise error). The third term is the error in density resulting when the smoothing 
function does not adequately represent the true path of the sphere. In this case, a bias 
error will result in density since the smoothing function will not fit the real perturbation 

in the path. 


NOISE ERROR TERMS 

The noise errors in vertical velocity (erf) and acceleration (c r z ) depend upon 

The noise present in the radar coordinates (a z ) 

The type of smoothing technique used 

The number of data points (N) used in the smoothing process 
The time spacing between consecutive data points (At) 

For an FPS-16 radar, cr z varies between 10 and 15 meters depending upon slant range 
and At is generally fixed at 0.1 second. 

ESTIMATION OF NOISE ERROR 

There are two methods of evaluating the noise error terms: 

(a) Consider an actual flight of a passive sphere tracked by two identical FPS-16 
radars. For an N and a smoothing function, density can be calculated for each of the two 
radar tracks. By calculating the RMS difference between the densities measured by the 
two radars, the noise error terms can be determined. Since the same bias error will 
appear in the density computations from each of the two radar tracks, differencing the 
densities determines only the noise error terms. 

(b) The noise error terms can also be calculated by formulas which directly relate 
and ag to N, At, a z , and the smoothing function. The formulas for polynomial 1 

smoothing functions of degrees one and two are given in the following equations: 

1 Polynomials were chosen as the proper class of smoothing functions. This deci- 
sion was based on previous work showing the polynomial yielding less noise error than 
other functions. 
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Noise Error Equations 

Velocity: \ 


Linear fit 

a 2 =• 

12 

A 



X 1 

n(n 2 - l) 

At 2 


Cubic fit 

a 2 = 

12 

( 7(3N 2 - 7) 2 (N - 4): 

4 


x 3 

_n(n 2 - l) 

(n+3): 

At 2 


Quadratic fit 

Acceleration: 
Linear-linear fit 

Cubic -linear fit 

Cubic -cubic fit 

Linear -cubic fit 


Same as linear fit o? = o? 

x 2 x l 


a? 


a 2 . 

X 11 M(M 2 - l) At 2 


12 _^3 
X 31 M(M 2 - l) At 2 


> ( 6 ) 


o2 = 

x 33 


c 13 


12 7(3M 2 - 7) 2 (M - 4)1 


M(M 2 - l) (M + 3)! 


a? 

x 3 


At 


12 7(3M 2 - 7) Z (M - 4).' 

m(m 2 - l) + (M + 3): 


1 


n 4 

X 1 


At: 


Quadratic (second derivative) a? = 

x 2 


720 

a x 


N 5 - 5N 3 + 4N 

At* 

J 


(Linear polynomial smoothing is defined as fitting a linear polynomial over N data points 
and assigning the slope of the fit to be the velocity at the midpoint ^^±-1 of the interval. 
Linear-linear smoothing to obtain acceleration is described as fitting N position points 
to a polynomial to obtain velocities and obtain acceleration from velocities in a like man- 
ner. A cubic -linear fit is described as fitting N position points to a cubic polynomial 
taking the slope at the midpoint as the velocity and fitting M of these velocities by a 
linear polynomial to obtain acceleration; similarly, for cubic-cubic and linear-cubic 
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smoothing techniques. Quadratic smoothing is defined as fitting a second-degree poly- 
nominal to position points and evaluating the first and second derivatives of the poly- 
nomial, at the midpoint, as the velocity and acceleration, respectively.) The validity 
of these formulas has been established by comparisons to RMS errors obtained by 

method (a). 


BIAS ERROR TERMS 

The bias errors in velocity and acceleration (Az,Az) depend upon 
The type of smoothing technique used 

The number of data points (N) used in the smoothing process 
The time spacing between consecutive data points (At) 

The true position field, which itself is a function of the balloon apogee 

ESTIMATION OF BIAS ERRORS 

For a given apogee, bias errors can be determined by following the flow chart of 
figure 4. Given a drag table and balloon apogee, by assuming the sphere fell in the 1962 
Standard Atmosphere, the equations of motion can be integrated to obtain the theoretical 
path of the sphere. The vertical position z, velocity z, and acceleration z are deter- 
mined by the theoretical trajectory. One now treats the z position as a function of time 
as though it were the radar coordinates and applies the smoothing routine using N data 
points and the degree polynomial P. The smoothed z, z, and z coordinates differ 
from the theoretical z, z, and z coordinates only because of the bias error resulting 
from the smoothing technique. (No noise has been introduced into the data.) The 
smoothed coordinates are then substituted into the equations of motion using the same 
drag table, and density is computed. The only difference between this computed density 
and the original input density, that is, the 1962 Standard Atmosphere density, is due to the 
bias error induced by the smoothing function. The percent bias error in density is then 
plotted as the ratio of the computed density to the standard density. 

DETERMINATION OF OPTIMUM SMOOTHING TECHNIQUE 

The choice of the optimum polynomial and interval for use in the high-altitude pro- 
gram was derived by using the above techniques. Initially, an escape altitude of 125 km 
was chosen. For each type double smoothing (linear -linear, cubic-linear, linear-cubic, 
cubic-cubic) the total error in density, defined as the square root of the sum of the noise 
and bias errors squared, was computed for all possible combinations of N and M. The 
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noise error was calculated by method (b) (the formulas as given in equations (6)) using 
At = 1/2 second and a z equals 15 meters. The bias error was computed as described 
in figure 4. Figures 5, 6, and 7 are plots showing the percent bias error (the deviation 
of the density ratio from 1) and the la confidence bands of the noise error about the 
bias. None of the combinations of the degree polynomials and smoothing intervals met 
the design requirements of 5% density accuracy at 100 km. A compromise was neces- 
sary to either reduce the altitude requirements, say to 95 km, maintaining the 5% design 
accuracy requirements or slackening the accuracy requirements to say 9% while main- 
taining the 100 km altitude range. The former choice was made. Careful analysis of 
total error plots for all combinations of degree polynomial and N and M resulted in 
the choice of the 19-21 linear-cubic combination as optimum. 

Total error plots were generated in the same fashion using a quadratic polynomial 
fit and its first and second derivatives for velocity and acceleration. The best smoothing 
interval for using a quadratic was determined to be 31 data points (fig. 7). In comparing 
the optimum quadratic and the optimum linear-cubic smoothing techniques, it is easily 
seen that the 19-21 linear -cubic produced significantly better results in the altitude 
region from 70 to 100 km. The probable explanation for this is as follows. By fitting 
two different functions, one to get velocity and the other to get acceleration, it is possible 
to partially compensate for, say, a positive bias in density due to a velocity error by 
using a different degree polynomial or different interval to generate accelerations which 
will produce a negative density bias. This advantage is not present when using a single 
function for smoothing. 

OPTIMUM FITTING FUNCTIONS TO ACQUIRE DENSITIES 

After careful consideration of the results of all the above analysis it was decided 
that the 19-21 linear-cubic smoothing was the optimum fitting function to acquire densities. 

SUMMARY OF DENSITY ERRORS USING OPTIMUM REDUCTION 

The total percent error in density resulting from a computation using the high- 
altitude ROBIN program with optimum smoothing cannot be precisely determined because 
of the unknown accuracy of the drag table and the occurrence of unknown vertical winds. 
However, the other contributing terms to density errors have been accurately determined 
using the optimum smoothing interval. They are less than 1 % bias error from 30 to 70 km 
with a noise error of 3% and a bias error from 3% to 5% from 70 to 95 km with a noise 
error less than 3%. 
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WINDS 


The equation used for computing wind with the falling sphere method is given as 

gx v BP\ 


W x = x 


(4 - W z )(x + C x - g x - 
= i 


z + c z - gz 


Sz v BP 


( 7 ) 


m 


The variables in equation (7) which may contribute significantly to wind errors are 
the horizontal and vertical components of velocity and acceleration. 

WIND ERROR EQUATIONS 

Having retained only those variables which significantly influence the error in a 
computed wind, the wind equation simplifies to 


W x = x - rr— 


z - g 

The first order of approximation to an error in W x resulting from the errors in 
the parameters is given by the following equation: 


( 8 ) 


dW x = dx — dx — dz + 


xz 


z - g 


z - g 


(z - g)' 


dz 


( 9 ) 


If the differential error components are considered as noise error with normal 
distribution, then by taking the variance of equation (9) the noise error in a wind calcula- 
tion is given as 


o o 
°W X = °k + 


z 

2 2 
(T.. 4- 

X 

2 0?+ 

xz 

z - g 

°x + 

Z - g 

°z + 

.(z - g)^ 


(10) 


where cr^, cr», o^., and a- are the noise errors in velocity and acceleration due to the 
noise in the radar data. To determine the bias of a wind measurement resulting from the 
bias (oversmoothed) velocity and acceleration measurements, equation (9) is again applied. 
Considering the component differentials as bias error, the square of the bias wind error 
is given as 

-.2 


a2w x = 


Ax - — Ax - — Az + - Az 


z - g 


- g (z - g)2 


( 11 ) 
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where the Ax, etc., refers to the bias error in the x velocity component. The total 

wind error ( 0 ™ \ is defined as the square root of the sum of the noise errror vari- 

V x total/ 

ance plus the bias square error and is given as 

^total = + i2Wx <12) 

The problem simply stated is to determine the type smoothing (degree) and smoothing 
intervals which minimize equation (12). As in the case of the density smoothing the noise 
error will decrease as the smoothing interval increases, and the bias error increases as 
the smoothing interval increases so that a minimum does exist for equation (12). 

MINIMIZATION OF TOTAL WIND ERROR 


The calculation of noise error for the wind computation uses equations (6). How- 
ever, the calculation of the bias error cannot use figure 4 because, with such a variety 
of wind profiles in nature, choosing one profile to be representative would not be realistic. 
Besides it was felt that the use of the following bias equations would be more precise: 


Bias Error Equations 

Position assumed 4th degree over N data points, i.e., ' 

x = Aq + A^ + A 2 t 2 + A 3 t 3 + Ajt 4 


Velocity assumed cubic over M data points, i.e., 

x = B 0 + Bit + B 2 t 2 + B 3 t 3 


Velocity: 
Linear fit 

Cubic fit 


Axi = 


A 3 At 2 (3N 2 - 7) 
20 


Ax 3 = 0 




Acceleration: 
Linear -linear 

Cubic -linear 
Cubic -cubic 


Ax ll =x 33 " hi 

BgAtj(3M 2 - 7) 


Ax 3 i = 


Ax 33 = 0 


20 


J 


(13) 
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To use these equations requires a knowledge of the position field of the balloon. This 
was accomplished by a separate program which utilizes experimental data to compute 
the coefficients A3 and A5. It is beyond the scope of this paper to describe that 
program. 


DETERMINATION OF OPTIMUM SMOOTHING FUNCTIONS 

The optimum double smoothing technique is that combination of degree polynomials 
(cubic-cubic, cubic-linear, linear-linear) and smoothing intervals (N and M) which 
gives the minimum total wind error. Plots of the total wind error for each type double 
smoothing and for N-M values of 51-43, 53-11, and 53-25 are presented in figures 8 to 
10. These are merely three illustrations of all the possible combinations for feasible 
values of N and M. After analyzing plots of the types illustrated by figures 8 to 10, it 
was determined that the 51-43 cubic -cubic smoothing provides optimum wind reduction. 

SUMMARY OF WIND ERRORS USING OPTIMUM REDUCTION 

With the type smoothing described above, the total wind error remains less than 
10 m/sec to altitudes of nearly 100 km. The amount of detail that can be observed at the 
very high altitudes is, however, limited because of the large altitude layer used in the 
smoothing. The frequency response curves presented as figure 11 indicate the detail 
that can be observed. Plotted in this figure is the ratio of the amplitude of a sinusoidal 
wave after passing through the smoothing filter to the true amplitude of the original. 

Each curve gives the ratio as a function of wavelength at a specific altitude. For example, 
at 90 km altitude the amplitude of a 10 km vertical wave would appear to be only 1/5 the 
amplitude in the reduced data. A 20 km vertical wave would retain 65 percent of its 
amplitude in the reduced data. As seen from figure 11, for 70 km, wavelengths less than 
10 km are largely destroyed so that only a mean wind profile can be ascertained. Below 
70 km, wavelengths of 5 km and less will appear in the reduced data. 

OPTIMUM PROGRAM DENSITY RESULTS 

A series of three Viper -Dart flights were flown at Eglin Air Force Base, Florida, 
on February 18, 1968, at 18:00, 19:00, and 20:00 zulu. Each flight was tracked by two 
FPS-16 radars. The flights are identified as Viper-Dart 11, 12, and 13 and the radars as 
radar 23 and 27. Figure 12 shows the density ratio (compared with the 1962 Standard 
Atmosphere) for each track of the three flights. For each of the flights, there is excellent 
agreement between the two FPS-16 radar tracks. The small differences in densities that 
are observed are commensurate with the noise errors predicted for the 19-21 linear-cubic 
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smoothing. There are, however, variations in density as observed from flights 1 hour 
apart, particularly in the altitude region from 62 to 54 km. The cause of these differ- 
ences is not known. Possible causes are 

(a) An actual time fluctuation in density 

(b) Special variation in density between the paths of the three spheres 

(c) Inaccuracy in the drag table being experienced at different altitudes for the 

three flights 

(d) A change in the vertical motions of the atmosphere 

These discrepancies in density, 1 hour apart, are not due to the inability of the 
radar to accurately track the spheres. Comparison of densities from the two tracks of 
the same balloon clearly rules this out. Nor are the density discrepancies thought to be a 
result of balloon collapse or elongation. All balloon collapse checks indicate the balloon 
is still spherically inflated to at least a 45 km altitude. 

OPTIMUM PROGRAM WIND RESULTS 

Figures 13 and 14 are plots of the W x and Wy components for each of the three 
flights. Both the W x and Wy components obtained from both tracks of Viper -Dart 11 
show nearly identical agreement. Viper -Dart 12 shows good agreement at altitudes below 
85 km but only fair agreement above. Viper -Dart 13 gives good agreement only to 84 km. 
The cause of this decrease in agreement which is beyond what one should anticipate from 
the total error plots for 51-43 cubic-cubic smoothing (fig. 8) has been investigated and 
the following results obtained. 

The tracks of radar 23 from both Viper -Dart 12 and Viper -Dart 13 flights show 
large oscillations which did not appear in the tracks of radar 27 from the same flights. 
Previous experience with FPS-16 tracking data indicates that the oscillations are prob- 
ably due to a low servo -bandwidth setting. The fact that radar 27 produces a smooth 
nonoscillating track indicates the oscillations are not real. Further investigation of the 
effect of the servo on tracking of passive spheres is in order. 

SUMMARY OF RESULTS 

Essentially there are three ways of determining the accuracy of the density and wind 
data: equations (6) and (13) and the model simulation outlined in figure 4. The assump- 
tions made in applying equations (6) and (13) and the model given in figure 4 are not pre- 
cisely met by the data but are exact enough for their purpose, which was to obtain the 
optimum smoothing technique. The use of double track flights to obtain errors only gives 
the noise error of the system. The dual tracking wind data shows that the radars need to 
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be well tuned in order to acquire winds, and that the tracking problem is more critical to 
winds than to density. Flights close together in time and space should yield system error 
information provided the time and space separation is small enough to rule out actual 
changes. One hour apart for density should be small enough to pick up the system errors. 
Comparison of the curves in figure 12 indicates that there are system errors in density. 
These system errors can conceivably come from the tracking radar, but they are larger 
than the equations predict. The system errors are believed to come probably from other 
elements of the system, sphere shape, or drag table. To strengthen this conclusion, the 
wind plots (figs. 13 and 14) have to be examined. The accuracy of the wind data does not 
depend upon the drag table or sphere shape (the last is almost true) but depends rather 
heavily upon the tracking radar. Figures 13 and 14 show that the wind repeatability is 
good, which proves that the tracking is good. Therefore, the conclusion is made that the 
discrepancy in density is not due to radar tracking but to other elements in the system. 
The dual tracked data agrees well with the noise errors predicted by equation (5) for den- 
sity and by equation (10) for winds. There is no other proof given here that the predicted 
bias errors are correct. This proof could come from simultaneous flights of different 
type sensors. 


TEMPERATURE AND PRESSURE 

Since the ROBIN is a density and wind sensor, the program optimized these vari- 
ables. Temperature and pressure errors fall where they may. The results of these 
parameters for the Viper -Dart flights 11, 12, and 13 are given in figures 15 and 16 without 
comment. 
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APPENDIX 


SYMBOLS 


Units are not given in the symbol list but any consistent set of units may be used in 
the equations. 

A cross-sectional area of sphere 

AMP amplitude 

A n ,B n coefficients in bias error equations, n = 0,1,2,... (see eqs. (13)) 

Cjj drag coefficient 

C X ,C Z Coriolis acceleration in x- and z-direction, respectively 

g gravitational acceleration 

g g gravitational acceleration in x- and z-direction, respectively 


K 

M 


m 

N 

t 

At 



v 


balloon constant 

number of data points used in acceleration smoothing process 
mass of sphere 

number of data points used in velocity smoothing process 
time 

time spacing between consecutive position data points 
time spacing between consecutive volume data points 
volume of balloon 

relative velocity of balloon, J^x - W x )^ + ^y - Wy^ + (z - W z 



1/2 
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APPENDIX 


W x ,W y 


W„ 


x,y,z 


z th ,z emp 


a 


Aq 

P 

p 0 


wind velocity in x- and y-direction, respectively 
vertical wind 

position coordinates of radar 
balloon velocity in x-direction 
balloon acceleration in x-direction 
balloon velocity in y-direction 
vertical velocity 
vertical acceleration 

theoretical and empirical vertical acceleration, respectively 

density gradient constant 

bias error of parameter q 

density 

initial density 

variance of parameter q 
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Drag coefficient as a function of Reynolds 
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Figure 5.- Ratio of smoothed to standard density (middle curve) using 19-21 linear-cubic 

smoothing and la noise error about bias. 
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smoothing and la noise error about bias. 




Figure 7.- Ratio of smoothed to standard density (middle curve) using 31 quadratic smoothing 

and la noise error about bias. 




Figure 8.- Total wind error plots for 51-43 cubic -cubic, cubic -linear, and 

linear -linear smoothing. 



Figure 9.- Total wind error plots for 53-11 cubic-cubic, cubic-linear, and 

linear -linear smoothing. 




Figure 10.- Total wind error plots for 53-25 cubic -cubic, cubic-linear, and 

linear -linear smoothing. 
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Figure 11.- Ratio of the amplitude of a sinusoidal wave as it would appear in the 
reduced data to its true amplitude at various altitudes. 51-43 cubic-cubic 
smoothing; escape altitude, 125 km. 





DENSITY RATIO 


VIPER-DART II 



Figure 12.- Ratio of density from Viper-Dart flights 11, 12, and 13 for 
each radar compared with 1962 Standard Atmosphere density. 

19-21 linear-cubic smoothing. 


VIPER-DART II 



Figure 13.- x-wind component from Viper -Dart flights 11, 12, and 13 for each radar 

51-43 cubic-cubic smoothing. 
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51-43 cubic-cubic smoothing. 




TEMPERATURE, KELVIN 



Figure 15.- Temperature from Viper-Dart flights 11, 12, and 13 for each radar 
compared with 1962 Standard Atmosphere temperature (dotted line) using 
19-21 linear-cubic smoothing. 
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Figure 16.- Ratio of pressure from Viper-Dart flights 11, 12, and 13 for 
each radar compared with 1962 Standard Atmosphere pressure. 

19-21 linear-cubic smoothing. 
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